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Liquid phase NMR is a general purpose test-bed for developing methods of coherent control 
relevant to quantum information processing. Here we extend these studies to the coherent control 
of logical qubits and in particular to the unitary gates necessary to create entanglement between 
logical qubits. We report an experimental implementation of a conditional logical gate between two 
logical qubits that are each in decoherence free subspaces that protect the quantum information 
from fully correlated dephasing. 



I. INTRODUCTION 




Experimental implementations of quantum information processors have reached sufficient complexity that it is now 
possible to experimentally explore the avoidance and correction of quantum errors by encoding quantum information. 
These information encoding methods include active techniques, like quantum error correcting codes (QECC)P, 0, 
an d dynamical decoupling 0, 0JE 0, and passive techniques, such as decoherence-free subspaces (DFS) 
noiseless subsystems (NS) [Rjj and topological schemes [l4|- To date the ability to store information 
Il7l ITsL llPJ ] , to perform universal quantum operations pfj , and to implement simple two qubit algorithms 
within a variety of logical encodings [2TL l22j| have been demonstrated. 

Entanglement, a uniquely quantum resource, enables many of the speed-ups afforded by QIP including many- 
body physics simulations [23, 13 exponential algorithmic enhancements [2j| , metrology j2|| and communication IpTl . 
Creation of entangled quantum states continues to drive experimental research in quantum information psl l29l l30l 
l3lj and has served as a benchmark for coherent control. Here we combine the two concepts of logical qubits and 
entanglement creation to prepare a pseudo-pure version of a Bell state between logical qubits. 

Control of encoded qubits must naturally respect the symmetries involved in the encoding. In the simplest case, this 
is achieved by having the control Hamiltonians commute with the noise generators 0] • When this is not experimentally 
possible or inconvenient, high fidelity control is achievable via modulation schemes that limit the encoded information's 
excursion out of protected subsystems to times short compared to the noise correlation time |32| . 

This report focuses on experimental implementations of the modulation sequences studied in |32^ for creating 
entanglement among logical qubits, specifically creating a logical Bell state between a pair of DFS qubits immune 
to collective dephasing. Furthermore, we perform our entanglement creation gate on two distinct initial states: (i) a 
pseudo-pure state effectively pure over the entire four qubit Hilbert space and (ii) a subsystem pseudo-pure state [see 
previous paper] . We also give an analysis of the quantum gate fidelities given these two input states and identify the 
largest errors in implementing this gate. 

II. LOGICAL BASIS ENCODING 

Using the open quantum system approach, we model the total Hamiltonian of our system and environment as : 

H = H s <8> t E + l s <g> He + Use (1) 

where Tis is the nuclear spin system Hamiltonian, He is the environment Hamiltonian and Hse describes the system- 
environment coupling. For this example we choose an encoding for a simple noise model: collective a z noise which 
corresponds to random fluctuations of the local magnetic field, B z . Defining the total angular momentum of the 
system as J z = X)i!=i °%i the interaction Hamiltonian is : 

H S e=iJz®B z . (2) 

The potential errors that the coupling to the magnetic fields can induce belong to the interaction algebra A z = 
{l, J z , J 2 , J^} [2(j. In the two spin case (N=2), the eigenspacc of the noise operator J z with eigenvalue is a 
C 2 x C 2 decoherence free subspace and can be used to encode one qubit of information. This DFS is thus spanned 
by the basis vectors 1 01} and |10). A natural encoding of a logical qubit l is given by: 

a\O) L +0\l) L ^a\Ol)+0\lO) (3) 



2 



The logical analogs of spin operators a x , a y . a z , and 1, which fully parametrize a single qubit, are: 



r-k J-\ 1* Z-Z. rr 1 ^ _2 S U 2. 



Furthermore, it shall be convenient to describe conditional logic using the logical idempotent operators: 

E£ & lL± er * = °lT<rj (5) 

In the four spin case (N=4) the J z eigenspace with eigenvalue is spanned by 6 basis vectors: |0011),|0101), 
|0110), 1 1001) , 1 1010) , and 1 1100) . Any four of these states can span a C 4 x C 4 subspace containing two logical qubits 
of information. We chose the basis |0101), |0110), 1 1001) , 1 1010) . In addition to being immune to fully correlated 
dephasing under J z , these states are also immune to pairwise collective dephasing under noise generators j 12 = <x* +CT 2 . 
and j]^ 4 = <7 Z + a 4 . The protected subspace is thus a tensor product space of two qubits of the form ©. 



III. IMPLEMENTING ENTANGLEMENT CREATION GATES ON A QUANTUM SYSTEM 

Given a fiducial state in the computational basis, applying a Hadamard gate and subsequently a CNOT gate, 
creates one of the four Bell states. The creation of a logical Bell state amounts to implementing logical versions of 
Hadamard and CNOT gates. The Hadamard gate on a logical qubit is specified, up to a global phase, as 

iU H = fZ+f* = e -if< e -*f <tf e if (6) 
v2 

The unitary operator for implementing a CNOT gate can be decomposed into a product of unitary operators of the 
form of single logical spin rotations and couplings of the 'ZZ' form: 

TT — T?~LL„2L i plL-a2L 

UCoNOT — a x + 1 

By expanding each of the exponentials above using the logical Pauli operators ijljl. many simplifications are possible. 
For instance, when a logical operator consists of a sum of commuting bilinear terms, it may suffice to drop all but one 
of the terms in the sum and add a constant scaling factor provided this simplified unitary has the same effect as the 
full unitary on a state within the logical encoding. One example is the isomorphism between a rotation of ^ about 
cr 4L and a ir rotation about cr^a 2 or CyC 2 ,: 



( cos(-)l - i sin(-)c7 4 a2) ( cos(-)l - i sin(-)o^) 



= cos 2 (-)l + sin 2 (-)a> 2 _ isin (0) 



(8) 



-ido 



= cos(0) — i sm(6)alcrl 
= cos(0)l - i sin((9)(cr+(7+ 
-- L cos(0)l - i sm(e)(a 1 + a 2 _ 



a\a 2 _ 



(9) 



cos(6>) - ism(0)<jl<j 2 
cos(6»)l -isin(0)( 
cos(0)l - isin(6»)(crV 2 



1 2 



a] a 2 



a 1 a 2 



a 1 a 2 



(10) 
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Here the last lines in II Oj l are all equivalent. The = L operation keeps only the terms that act within the logical 
subspace. Explicitly, the logical basis states are eigenstates of the a\a1 operator with eigenvalue -1; likewise, the 
eigenvalues of o\xr± are all zero. Using the following substitutions the net unitary of a Hadamard gate and a CNOT 
gate can be decomposed into four separate logical operations (up to a global, unobservable phase), each consisting of 
a ? rotation about a single bilinear term: 



UcoNOtUh & e^UiUaUaUi 



(11) 



Ui 

u 3 



■ 4 "v"v 
t Tt _2 3 



u 2 

h 



(12) 



The <7 • a interaction present in the internal Hamiltonian provides a means to perform rotations, but this is 
efficient only when the chemical shift differences (|cjj — Wq\) are suppressed. When the chemical shift differences are 
non- negligible a Carr-Purcell (CP) style sequence can be used to engineer a a ■ a interaction j33, |3i|. By applying 
collective ir rotations on a pair of spins comprising a logical bit, the chemical shift terms can be averaged to zero, 
retaining the ZZ term to first order. Staggered n rotations on the spins of the other logical bit are applied in a 
manner to refocus all operators in the Hamiltonian associated with these spins (see Fig. IIIH . 





f | A | A | A | A | 

90x x x x 90x 

I I I I I 

90x x x x 9 0x 



_L 



J_ 



zxz 



FIG. 1: From left to right: entangling circuit on logical qubits and corresponding logical pulses and the puulse sequence 
implementing the a^ L logical rotation on the physical qubits. 



To generate a a x coupling, we first generate the zeroth order average Hamiltonian Z\ 
employing a CPdike sequence: 



s-Wo* mm by 



A12 - 7T) 



-Al2-7T) 



-A 12 -7T)1 



Ai 2 -^ 



(13) 



(Note the notation 0)* = exp[— i| ^ fe cr^].) The effective unitary propagator to zeroth order is = e 4lA ^z 12 ^ 

By applying "collective" ^ rotations about x before the above sequence and a -| rotation about — x after, the 
average Hamiltonian is transformed into &y&y which has the same action as a <j\ l . Similarly, if the rotation axes are 
separated by ^ (i.e. an y phased pulse on one of the spins in the logical pair and a x phased pulse on the other), 
an operator isomorphic to a a y L is achieved. The logical "two-body" interaction - a\ L cr^, - acting on the encoded 
subspaccs is isomorphic to o\a\. We obtain this by using l|13|) and replacing spin 2 with spin 3. Each of the four 
rotations in (|12|l can be generated in this manneryielding an overall sequence: 



1,2 .^1,2,3,4 7T\ 3 ' 4 7T\ 3 7T\ 4 

* - Zl2 ~2) x " Z34 ~2J 2 " 223 ~ 2 J, 2 J, 



Z.34- 



(14) 



where Zjk denotes the subsequence generating Zjk and setting A^fc = (2 J. 



IV. EXPERIMENT 



Our quantum system is comprised of the four 13 C spins (I = \) of isotopically labeled crotonic acid (also known 
as trans-3-butenoic acid) |3(i| |. This molecule contains nine magnetically active nuclei in total; the total spin system 
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Hamiltonian takes the form: 

Hinternal = £ ~ ^Vz + T^Y.^ ~ W oV* + Y g ' & + (15) 

ieC fceH i,j£.C;i<j 

fe,iGH;fc<; j,k;j<k 

where u>q (uq~) is the rotating frame frequency near the 13 C ( 1 H) Larmor frequency, the u)i are the chemical shifts 
of the 4 carbon nuclear spins, the uik are the chemical shifts of the 5 hydrogen nuclear spins, and the 6 Jy (10 
Jki) are the scalar coupling constants between two carbon (two hydrogen) spins (as usual, H = 1). We are mainly 
concerned with coherently controlling the carbon subsystem of spins and seek to suppress the proton subsystem. 
As the heteronuclear scalar couplings (terms Jj t k) are the only means of mixing the two subsystems, a broadband 
decoupling sequence modulating the proton spin system effectively removes this coupling during the experiment. In 
practice, decoupling the proton spin system is equivalent to saturating the populations of the proton spins. One 
potential artifact of this approach is the introduction of transient nuclear Overhauser effects (NOE) [371 HE Hil ■ 



A. Pseudo-pure states over the entire Hilbert Space 



Starting with the equilibrium density matrix of the four spin system, p « 1/iV — ^Q2±_ cr i)j we prepared our system 
in the dual DFS ground state, |0101) = |00)l using spatial averaging techniques |4^,|4^. As shown in Figure |2| the 
preparation of the encoded logical state is complex and requires a preparation time of 0.1186 seconds (s) relative to 
the total experiment length of 0.1662 s The T2's of the carbons are all greater than 500 milliseconds (ms) 01 and 
therefore spin-spin relaxation is unimportant over the length of the experiment. Since pseudo-pure state preparation is 
a non-unitary, completely positive map, a loss of observable signal is expected. In our implementation of pseudo-pure 
state preparation the signal is roughly 2/13 that of the equilibrium state. 



c 4 _ 



71/2 



71/2— * 



it : 



n/2 



FIG. 2: Circuit for the preparation of the pseudo-pure state. We represent single qubit rotations by square boxes, controlled 
rotations by closed circles on the controlling qubit linked to the applied rotation on the controlled qubit; swaps gate by two 
crosses on the swapped qubits, connected by a vertical line; non unitary operations (gradients) by double vertical bars. Notice 
the number of controlled operations, each requiring a time of the order of the inverse coupling strength, and SWAP gates, each 
requiring three times the inverse coupling strength. The single-qubit rotations above have the values a m t^tv, (3 ~ |-7T, and 
7 ~ j it which account for the scaling of the signal-to- noise compared to the equilibrium state. 

All qubit rotations (selective, semi-selective, or collective to all spins) were created using robust strongly-modulating 
pulses (SMP) [44ll45| by maximizing the gate fidelity of the ideal propagator to the simulated propagator. Furthermore, 
as the radio frequency (RF) control fields are inhomogeneous over the sample, we maximize the effective gate fidelity, 
averging over a weighted distribution of RF field strengths. Pulse lengths for this system range from 200-800 ps; 
the simulated gate fidelities for any individual pulse are greater than 0.99%. The dominant source of residual errors 
in a typical SMP come from two-body terms of the form c^cr^. After many pulses these small residual errors 
accumulate, but the net effect can be partially suppressed by adjusting the delay time between pulses |46l | to optimize 
the overall gate fidelity or state correlation. The required ^ rotations about the logical operator axes were obtained 
by implementing the sequence in l|13|) and setting Ajk = 2j~Z' ^ n principle, we would like to repeat the sequence N 

times and scale A by a factor of in order to induce a rapid refocusing of the noise; however, in practice A is limited 
by the length of the semi-selective pulses. For all of the logical rotations implemented experimentally, we use N=l. 

The density matrices of the spin system after the initial pseudo-pure preparation and after the entanglement creation 
were reconstructed using state tomography 0, which involves applying 18 readout pulses to obtain coefficients for 
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the 256 operators comprising a complete operator basis of the four spin-i 13 C nuclei. 

B. Subsystem pseudo-pure states 

As shown in the accompanying paper |4S[ . for mixed-state ensemble quantum information processing there are 
advantages to requiring the initial state to be pure only over the subsystem containing the relevant quantum informa- 
tion. Here we implement the entangling operation over logical qubits using the double DFS initial state described in 
|48j . This state (Figure 0J can be prepared in half the time of the full pseudo-pure state (0.0568s) and has a smaller 
loss of signal (2/3 compared to 2/13). 




FIG. 3: Circuit for the preparation of the subsystem pseudo-pure state. We use the same convention as in Figure [5] 



C. Density Matrix Reconstruction 

The reconstructed density matrices are shown in Figures 01 and [S] where the vertical axis shows the normalized 
amplitude and the horizontal axes label the basis states in the computational basis (i.e. 1 0000) , 1 0001) ,...). The effects 
of decoherence can be qualitatively seen in the final state as an attenuation of the off-diagonal terms of the Bell state: 

100)^-111) j. 




(a) Subsystem Pseudo-Pure State 



(b) Logical Bell State 



FIG. 4: Density matrices for the initial pseudo-pure state over only the logical subspace (a) and the corresponding Bell-State 
(b). The darker part indicates the states in the logical subspace. In the case of the subsystem pseudo-pure stat es, the division 
of the logical subspace allows for the other areas of Hilbert space to be mixed. 



(a) Pseudo-Pure 



(b) Logical Bell State 



FIG. 5: Density matrices for the initial pseudo-pure state over the entire Hilbert space (a) and the corresponding logical Bell 
State (b). The basis states comprising the dual DFS have been darkened. 



V. DISCUSSION 
A. Metrics of Control 

The quantum process associated with the encoded entangling operation can be specified by the general map: 

£(p) = 5>np4 (16) 

where A p are Kraus operators describing the experimental implementation of the encoded operation, and p is defined 
over the entire Hilbert space. If the process contains no error only a single Kraus operator defines the operation: 
Aq = U. The correlation of two quantum states is defined as: 

c = Tr{p th Pexp} , 17 x 
VTr{p 2 th }Tr{pjJ 

which defines the relative closeness of states in Hilbert space, taking into account attenuation due to decoherent or 
incoherent processes. Here pth = U Pi n U' , p exp = £(pi n ), and pi n define the theoretical, experimental and input states 
respectively. The closeness of the theoretical state to the experimental state measures how well we have implemented 
the logical entangling operation for a particular input state. We reserve the symbol C for designating the correlation 
that does not take into account the attenuation of the input state; in this measure pi n is replaced by p exp . 

When the relevant quantum information resides only in a subspace, we can also measure the correlation inside this 
subspace by applying the projectors P\_ spanning the subspace to the state [igj : 

sjTr{(P L p„P L y}Tr{{P L p, h P L )*} 

This measure is important when creating subsystem pseudo-pure states, where nonzero contributions to the density 
operator outside the logical subspace can degrade C. If the input state of this process is a full pseudo-pure state 
(FPPS), then pi n = P\_pi n P\_- Furthermore, the result of a logical operation acting on a FPPS adheres to the relation: 
pth = PiPthPi- Under these conditions the correlation C|_i_ with the logical ideal state is just the total correlation C. 

Finally, we can consider the correlation of the state after (i) an ideal decoding operation, Uo cc which transforms 
the subspace structure of the Hilbert space back to a tensor product structure and (ii) a partial trace over the ancilla. 
This examines how the entanglement created under two logical qubits transfers back to physical qubits. For the 
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Quantum state 


Csim 


f" 

Exp 


C*Exp 


fil f 1 fit f i 

\ 1 wLL ^rw ^Dcc 
i_ i_ i ) <_. u ^ ^ 


Full pseudo-pure (FPPS) 

Full pseudo-pure Bell state (FPPBS) 


0.950 
0.747 


0.685 


0.884 
0.531 


0.981 0.986 
0.904 0.590 0.852 0.641 


Subsystem pscudo-purc (SPPS) 
Subsystem pseudo-pure Bell state (SPPBS) 


0.984 
0.919 


0.836 


0.965 
0.725 


0.988 0.992 
0.849 0.608 0.974 0.823 



TABLE I: Experimental and simulated data for the implementation of encoded Bell state propagator. Experimental errors 
of ~ 4% can be attributed to systematic errors in the fitting algorithm used to reconstruct the density matrix from NMR 
spectral data. The correlations of the ideal logical Bell state are reported with respect to the following states: a simulation of 
the spin system with distribution of RF powers (Csim), the experimental tomography reconstruction of the state in the entire 
Hilbert space (Ce xp ), the the experimental reconstruction projected onto the logical subspace (Cll). Cd cc compares the ideal 
state and the reconstructed state assuming ideal decoding out of the logical subsystem and tracing over the ancillae. Both the 
figures reported for the Full pseudo-pure and Subsystem pseduo-pure Bell states reflect attenuation due to relaxation processes. 
Primes (/) indicate unattenuated correlations. 

encoding chosen in this paper U^ cc is: 

Ubcc = ClNOT 2 ■ C^NOT 4 

= (Elal+EH 2 )(Eia$ + E 3 _t*) (19) 

C Dcc has the same form as ([T7| with p th , p exp , and p in replaced by pf h , p d exp , and pf n where = tr a {U Dcc p[ ]U£, ec } 
and a is the subsystem spanned by the ancilla qubits. These values arc displayed in Table [I] 

Since the same pulse sequence implementing the entanglement creation was used for both initial states, differences 
in the correlation are due the unwanted evolution of initial errors in state preparation. These differences reveal 
several key features of the experimental implementation. In comparing the correlation of the simulated resutls to the 
experiment, we see the smallest deviations in the case of the subsystem pseudo-pure state (SPPS), attributing this to 
the benefits of a shorter preparation sequence. We also note that T\ and T 2 processes, the NOE and RF Transients 
are not contained in the simulation model. It follows that the longer and more complex preparation sequence for the 
full pseudo-pure state (FPPS) should result in a larger deviation of the correlation. For the simulation of the full 
pseudo-pure Bell state (FPPBS) and the subsystem pseudo-pure Bell state (SPPBS), a reduction from unity is due 
solely to incoherent (RF inhomogenities) and coherent errors; the experimental correlations are further affected by 
the processes stated above, as well as T 2 processes. 

In comparing the experimental logical subspace correlation (Cll) to C, we can better understand the source of 
error. For instance, for the SPPS, FPPS, and FPPBS, Cll is greater than C exp . This difference represents errors in 
C exp due to unwanted, nonzero contributions to the density operator outside the logical subspace since C ~ ^ k a^C^ 
where Cjk represent correlations for blocks of the density operator and a is a weighting factor |48j . However, in the 
case of the SPPBS, Cll is less than C exp . Here o^l = 0.40, «rr = 0.58, and Crr = 0.822. Thus by removing the 
Crr contribution from C exp we get a more complete view of the control over the logical qubits. Consequently for 
this state we controlled the information in the non-logical subspace better than the information inside the subspace 
of interest. We believe the source of this difference to be T 2 decoherence, as the majority of information in the 1Z 
subspace exists as population terms (er 2 , er 2 er z , etc). In the case of fully correlated noise, the characteristic decoherence 
time for any state with J z = 0, T 2 , should approach the characteristic relaxation time of the population terms (Ti). 
For this implementation, the T\ times for any individual spin are greater than 2 seconds, whereas the shortest T 2 is 
1.1 seconds which is consistent with our state tomography measurements. 

Lastly, we comment on the differences between the Cll and Cdoc- Both of these metrics assess the closeness of two 
states in a C 2 (£> C 2 space, but in a different manner. In the absence of relaxation, the decoded correlation assesses how 
well the implemented entangling operation respects the subspace structure of the system, namely U en t = U en t L ®U en t R - 
Under the decoding, this operation should only create entanglement among two physical qubits and thus the ancillae 
can be traced over. To the extent that U en t does not respect the subspace structure, the decoding will in general 
introduce entanglement between the informational qubits and the ancillae, thus resulting in a loss of purity after a 
partial trace. In the case of the FPPBS, Cd c is larger than Cll due to a loss of purity of the state; however, Cq 6C is 
smaller than C LL . This shows that the sensitivity to coherent unitary errors and incoherent/decoherent errors of the 
two metrics differs slightly. In the case of the SPPBS, the correlation and unattenuated correlation of the decoded 
metric are greater than their respect logical projection metrics. This can be attributed to the contributions to p' exp 
from state preparation. (Recall that for a SPPS, the decoded state is not pure on the information carrying qubits, even 
in the ideal case.) These contributions are less sensitive to errors during the quantum operation and thus augment 
the value of the decoded correlation. 
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Quantum State tr{p 2 exp ) tr(p 2 slm ) 

Pseudo-pure 0.781 0.812 

Logical Bell State 0.470 0.556 

TABLE II: The above values correspond to the purity of the initial and final states for the experimentally derived density 
operator and that obtained by simulating the four carbon system with RF inhomogeneity. Note that the purity of the initial 
state is not unity, as small errors to the eigenspectrum of p are quadratically amplified by the trace. The loss of purity from 
the initial to the final state for the simulated gates does not take into account decoherence; thus, any loss of purity must arise 
from the incoherence due to RF inhomogeneity. 



VI. CONCLUSION 

Using a two-physical qubit which protects against collective dephasing, we have shown how to implement quantum 
gates between logical qubits using effective Hamiltonians. It is important to stress that when selecting a protection 
scheme against decoherence the ability to encode a physical qubit into a logical qubit with high fidelity is not sufficient 
for computation. The structure of the external and natural Hamiltonian plays an important role in the control of 
logical qubits, as the operators needed to implement gates may not be present and generating them may drive 
the information out of the subsystem. For large systems with significant symmetry (like quantum dots under the 
exchange interaction 0) or exceedingly small systems [2£j, the structure of the natural Hamiltonian can provide the 
logical operations in itself. However, for systems of intermediate size (most relevant to the present implementations of 
quantum information processors) implementing quantum gates among logical qubits requires both a precise knowledge 
of the system Hamiltonian and a set of control parameters large enough to ensure no leakage from the protected 
subsystem or subspace. For example, if our four-qubit system were composed of two protons and two carbons, each 
individual species could be modulated separately, thus doubling the number of parameters in the control Hamiltonian 
and limiting the leakge of the information from the subspace. 

Finally, our selection of logical qubits comprised of only two physical qubits limits our logical operations for single 
qubit and two-qubit interactions to only "two-body" operators. If instead we were to attempt a repetition of the 
experiment where the logical qubits were encoded under the three qubit noiseless subsystem |l(j , the single qubit and 
two-qubit rotations would involve "three-body" operators |50l | - quite unlikely to be found in a natural Hamiltonian. 
In such a scenario, the ability to implement logical operations would necessarily need to come from a modulation 
sequence, appropriately chosen to avoid leakage from the subsystem. 

VII. ACKNOWLEDGEMENTS 

This work was supported in part by the National Security Agency (NSA) under Army Research Office (ARO) con- 
tract numbers W911NF-05-1-0469 and DAAD19-01-1-0519, by the Air Force Office of Scientific Research, and by the 
Quantum Technologies Group of the Cambridge-MIT Institute. We also thank N. Boulant for useful discussions. 



[1] P. W. Shor, Phys. Rev. A 52, R2493 (1995). 

[2] A. R. Calderbank and P. W. Shor, Phys. Rev. A 54, 1098 (1996). 

[3] D. Gottesman, Phys. Rev. A 54, 1862 (1996). 

[4] A. M. Steane, Proc. R. Soc. Lond. A. 452, 2551 (1996). 

[5] A. M. Steane, Phys. Rev. Lett. 77, 793 (1996). 

[6] L. Viola and S. Lloyd, Phys. Rev. A 58, 2733 (1998). 

[7] L. Viola, S. Lloyd, and E. Knill, Phys. Rev. Lett 83, 4888 (1999). 

[8] L. Viola, E. Knill, and S. Lloyd, Phys. Rev. Lett. 82, 2417 (1999). 

[9] L. Viola and E. Knill, Phys. Rev. A 66, 12307 (2002). 
[10] P. Zanardi and M. Rasetti, Phys. Rev. Lett. 79, 17 (1997). 
[11] L.-M. Duan and G.-C. Guo, Phys. Rev. Lett. 79, 1953 (1997). 
[12] D. A. Lidar, I. L. Chuang, and K. B. Whaley, Phys. Rev. Lett. 81, 2594 (1998). 
[13] E. Knill, R. Laflamme, and L. Viola, Phys. Rev. Lett. 84, 2525 (2000). 
[14] A. Y. Kitaev, Annals of Physics 303, 2 (2003). 

[15] P. G. Kwiat, A. J. Berglund, J. B. Altepeter, and A. G. White, Science 290, 498 (2000). 

[16] E. Fortunate, L. Viola, M. Pravia, E. Knill, R. Laflamme, T. Havel, and D. Cory, Phys. Rev. A 67, 62303 (2003). 



9 



[17] D. Kielpinski, V. Meyer, M. A. Rowe, C. A. Sackett, W. M. Itano, C. Monroe, and D. J. Wineland, Science 291, 1013 
(2001). 

[18] H. Haffner, F. Schmidt-Kaler, W. Hansel, C. F. Roos, T. Kober, M. Chwalla, M. Riebe, J. Benhelm, U. D. Rapol, C. Becher, 

et al., App. Phys. B 81, 151 (2005). 
[19] C. Langer, R. Ozeri, J. D. Jost, J. Chiaverini, B. DeMarco, A. Ben-Kish, R. B. Blakestad, J. Britton, D. B. Hume, W. M. 

Itano, et al., Phys. Rev. Lett. 95, 060502 (2005). 
[20] E. Fortunate, L. Viola, J. Hodges, G. Teklemarian, and D. Cory, New J. Phys 4, 5.1 (2002). 
[21] M. Mohseni, J. S. Lundeen, K. J. Resch, and A. M. Steinberg, Phys. Rev. Lett. 91, 187903 (2003). 
[22] J. Ollerenshaw, D. Lidar, and L. Kay, Phys. Rev. Lett. 91, 217904 (2003). 
[23] R. Feynman, Foundations of Physics 16, 507 (1986). 
[24] S. Lloyd, Science 273, 1073 (1996). 

[25] P. W. Shor, in Proceedings 35th Annual Symposium on Foundations of Computer Science (1994), pp. 124-134. 
[26] P. Cappellaro, J. Emerson, N. Boulant, C. Ramanathan, S. Lloyd, and D. G. Cory, Phys. Rev. Lett. 94, 020502 (2005). 
[27] C. H. Bennett and G. Brassard, in Proceedings of IEEE International Conference on Computers Systems and Signal Processing 
(1984), pp. 175-179. 

[28] C. Negrevergne, T. S. Mahesh, C. A. Ryan, M. Ditty, F. Cyr-Racine, W. Power, N. Boulant, T. Havel, D. G. Cory, and 

R. Laflamme, Phys. Rev. Lett. 96, 170501 (2006). 
[29] J.-S. Lee and A. K. Khitrin, Phys. Rev. A 70, 022330 (2004). 

[30] D. Leibfried, E. Knill, S. Seidelin, J. Britton, R. B. Blakestad, J. Chiaverini, D. B. Hume, W. M. Itano, J. D. Jost, 

C. Langer, et al., Nature 438, 639 (2005). 
[31] H. Haffner, W. Hansel, C. F. Roos, J. Benhelm, D. Chek-al kar, M. Chwalla, T. Korber, U. D. Rapol, M. Riebe, P. O. 

Schmidt, et al., Nature 438, 643 (2005). 
[32] P. Cappellaro, J. H. Hodges, T. F. Havel, and D. G. Cory, J. Chem Phys. 125, 044514 (2006). 
[33] H. Carr and E. Purcell, Phys. Rev 94, 630 (1954). 
[34] U. Haeberlen and J. S. Waugh, Phys. Rev. 175, 453 (1968). 

[35] U. Haeberlen, High Resolution NMR in Solids: Selective Averaging Academic Press Inc. (1976). 

[36] E. Knill, R. Laflamme, R. Martinez, and C.-H. Tseng, Nature 404, 368 (2000). 

[37] A. W. Overhauser, Phys. Rev. 91, 476 (1953). 

[38] A. W. Overhauser, Phys. Rev. 92, 411 (1953). 

[39] T. Carver and P. Slicher, Phys. Rev. 92, 212 (1953). 

[40] I. Solomon, Phys. Rev. 99, 559 (1955). 

[41] D. Cory, M. Price, and T. Havel, Physica D 120, 82 (1998). 

[42] D. G. Cory, A. F. Fahmy, and T. F. Havel, Proc. Natl. Acad. Sci. U.S.A. 94, 1634 (1997). 
[43] N. Boulant, T. F. Havel, M. A. Pravia, and D. G. Cory, Phys. Rev. A 67, 42322 (2003). 

[44] E. Fortunato, M. Pravia, N. Boulant, G. Teklemariam, T. Havel, and D. Cory, J. Chem. Phys. 116 (17), 7599 (2002). 
[45] M. A. Pravia, N. Boulant, J. Emerson, A. Farid, E. M. Fortunato, T. F. Havel, R. Martinez, and C. D. G., J. Chem Phys. 
119, 9993 (2003). 

[46] E. Knill, R. Laflamme, R. Martinez, and C. H. Tseng, Nature 404, 368 (2000). 

[47] I. L. Chuang, N. A. Gershenfeld, M. G. Kubinec, and D. Leung, Proc. R. Soc. Lond. A. 454, 447 (1998). 

[48] P. Cappellaro, J. S. Hodges, T. F. Havel, and D. G. Cory, quant-ph/0610269 (2006). 

[49] D. DiVincenzo, D. Bacon, J. Kempe, G. Burkard, and K. B. Whaley, Nature 408, 339 (2000). 

[50] C. H. Tseng, S. Somaroo, Y. Sharf, E. Knill, R. Laflamme, T. F. Havel, and D. G. Cory, Phys. Rev. A 61, 012302 (1999). 



